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Quantum Information Processing relies on coherent quantum dynamics for a precise control of 
its basic operations. A swapping gate in a two-spin system exchanges the degenerate states It, J.) 
and ||, T). In NMR, this is achieved turning on and off the spin-spin interaction b = AE that 
splits the energy levels and induces an oscillation with a natural frequency AE/h. Interaction 
of strength K/tse, with an environment of neighboring spins, degrades this oscillation within a 
decoherence time scale t^. While the experimental frequency uj and decoherence time were 
expected to be roughly proportional to b/h and tse respectively, we present here experiments that 
show drastic deviations in both ui and t^. By solving the many spin dynamics, we prove that the 
swapping regime is restricted to A_Etse > ft. Beyond a critical interaction with the environment 
the swapping freezes and the decoherence rate drops as 1/tj, oc (b/h) 2 tq E . The transition between 

quantum dynamical phases occurs when ui oc \J (b/h) 2 — (k/r SE ) 2 becomes imaginary, resembling 

an overdamped classical oscillator. Here, < k 2 < 1 depends only on the anisotropy of the system- 
environment interaction, being for isotropic and 1 for XY interactions. This critical onset of a 
phase dominated by the Quantum Zeno effect opens up new opportunities for controlling quantum 
dynamics. 



I. INTRODUCTION 

Experiments on quantum information processing^ in- 
volve atoms in optical traps^, superconducting circuits^ 
and nuclear spins^ among others. Typically, the system 
to be manipulated interacts with an environmenliSiSiLS 
that perturbs it, smoothly degrading its quantum dy- 
namics with a decoherence rate, proportional to 
the system-environment (SE) interaction H/tse- Strik- 
ingly, there are conditions where the decoherence rate can 
become perturbation independent £ This phenomenon is 
interpretedi2*ii*i£ as the onset of a Lyapunov phase, 
where 1/t0 = min[l/rsE,A] is controlled by the sys- 
tem's own complexity A. Describing such a transition, 
requires expressing the observables (outputs) in terms of 
the controlled parameters and interactions (inputs) be- 
yond the perturbation theory. We are going to show that 
this is also the case of the simple swapping gate, an essen- 
tial building block for quantum information processing, 
where puzzling experiments require a substantially im- 
proved description. While the swapping operation was 
recently addressed in the field of NMR in liquidsi 3 *^ 
with a focus on quantum computation, the pioneer exper- 
iments were performed in solid state NMRi£ by Miiller, 
Kumar, Baumann and Ernst (MKBE). They obtained a 
swapping frequency uj determined by a two-spin dipolar 
interaction 6, and a decoherence rate 1/ (2r^) = R that, 
in their model, was fixed by interactions with the envi- 
ronment 1/(2tse)- This dynamical description was ob- 
tained by solving a generalized Liouville-von Neumann 
equation. As usual, the degrees of freedom of the en- 
vironment were traced out to yield a quantum master 
equation^ More recent experiments, spanning the in- 
ternal interaction strength^ hinted that there is a crit- 
ical value of this interaction when a drastic change in 



the behavior of the swapping frequency and relaxation 
rates occurs. Since this is not predicted by the standard 
approximations in the quantum master equation^ this 
motivates to deepen into the physics of the phenomenon. 

In this paper, we present a set of 13 C- 1 H cross- 
polarization NMR data, swept over a wide range of a 
control parameter (the ratio between internal interac- 
tions and SE interaction strengths). These results clearly 
show that the transition between the two expected dy- 
namical regimes for the 13 C polarization, an oscillating 
regime and an over-damped regime, is not a smooth 
cross-over. Indeed, it has the characteristics of critical 
phenomena where a divergence of the oscillation period 
at a given critical strength of the control parameter is 
indicative of the nonanalyticity of this observable .iSiisi 
The data are interpreted by solving the swapping dy- 
namics between two coupled spins (qubits) interacting 
with a spin bath. With this purpose the environment is 
represented as a stroboscopic process. With certain prob- 
ability, this instantaneous interaction interrupts the sys- 
tem evolution through measurements and/or injections. 
This simple picture, emerging naturally22iSi from the 
quantum theory of irreversible processes in the Keldysh 
formalismjSSiS^ enables us to distinguish the interaction 
parts that lead to dissipation from those giving pure de- 
coherence. Within this picture, the overdamped regime 
arises because of the quantum Zeno effect ^iSiSii i.e. en- 
vironment "measures" the system so frequently that pre- 
vents evolution. The analytical solution confirms that 
there is a critical value of the control parameter where 
a bifurcation occurs. This is associated with the switch 
among dynamical regimes: the swapping phase and the 
Zeno phase. In consequence, we call this phenomenon a 
Quantum Dynamical Phase Transition. 
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II. EXPERIMENTAL EVIDENCE 

Our cross polarization experiments exploit the fact 
that in polycrystalline ferrocene Fe(C5H 5 ) 2 , one can se- 
lect a pair of interacting spins, i.e. a 13 C and its directly 
bonded X H, arising on a molecule with a particular ori- 
entation. This is because the cyclopentadienyl rings per- 
form fast thermal rotations (ss ps) around the five-fold 
symmetry axis, leading to a time averaged 13 C- 1 H in- 
teraction. The new dipolar constant depends^ only on 
the angle 9 between the molecular axis and the external 
magnetic field B and the angles between the internu- 
clear vectors and the rotating axis, which in this case are 
90°. Thus, the effective coupling constant is 



other protons that constitute the environment Their 
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where 7's are the gyromagnetic factors and r HC the in- 
ternuclear distance. Notice that b(9) cancels out at 
the magic angle 9 m ~ 54.74°. As the chemical shift 
anisotropy of 13 C is also averaged by the rotation and 
also depends on 9 as (3 cos 2 9 — I s ) , it is straightforward 
to assign each frequency in the 13 C spectrum to a dipo- 
lar coupling 6. Thus, all possible b values are present 
in a single polycrystalline spectrum. The swapping in- 
duced by b is turned on during the "contact time" t, 
when the sample is irradiated with two radio frequencies 
fulfilling the Hartmann-Hahn condition^ Experimental 
details have been given elsewhere^ At t = 0, there is 
no polarization at 13 C while the 1 H system is polarized. 
The polarization is transferred forth and back in the 13 C- 
1 H pairs while the other protons inject polarization into 
these pairs. In Fig. we show the raw experimen- 
tal data of 13 C polarization as a function of the contact 
time and b{9). The polarizations have been normalized 
to their respective values at the maximum contact time 
(3 ms) for each 9 when it saturates. It can be appreci- 
ated in the figure that the oscillation frequency is roughly 
proportional to |6| , showing that this is the dominant in- 
teraction in the dynamics. This is consistent with the 
fact that the next 13 C- 1 !! coupling strength with a non- 
directly bonded proton is roughly 6/8 and, as all the in- 
tramolecular interactions, also scales with the angular 
factor (3 cos 2 6- l). 

A noticeable feature in these experimental data is the 
presence of a "canyon", in the region |6| < 2 kHz, where 
oscillations (swapping) disappear. The white hyperbolic 
stripes in the contour plot at the bottom evidence a swap- 
ping period 2it/lo that diverges for a non-zero critical in- 
teraction. This divergence is the signature of a critical 
behavior. 

The standard procedure to characterize the cross 
polarization experiment in ferrocene and similar 
compounds^ is derived from the MKBE model. There 
the 13 C polarization exchanges with that of its directly 
bonded '■H, which, in turn, interacts isotropically with 
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where the decoherence rate becomes determined by the 
rate of interaction with the environment 1/(2tse) — ► 
1/ (2r^) = R, while the swapping frequency is given by 
the two-spin dipolar interaction, b/h — > lo. A dependence 
of the inputs b and tse on 9 should manifest in the observ- 
ables lo and T0. However, working on a polycrystal, each 
t~se(9) value involves a cone of orientations of neighbor- 
ing molecules and a rough description with single average 
value for the SE interaction rate is suitable. 

We have performed non-linear least square fittings of 

• • MKBE . . 

the experimental points to the equation P (t) tor the 
whole 13 C spectra of ferrocene in steps of w 80Hz and 
contact times ranging from 2/zs to 3ms. The 1/ (2r^) 
and lo parameters obtained from these fits are shown as 
dots in Fig. [21 The proportionality of the frequency 
with b for orientations that are far from the magic angle 
is verified. In this region a weak variation of l/(2r^) 
around 2.2kHz reflects the fact that 1/(2tse) does not 
depend on 9. A drawback of this simple characterization 
is that it tends to overestimate the width of the canyon 
because of limitations of the fitting procedure when Eq. 
© is used around the magic angle. 

In spite of the MKBE theoretical prediction, one ob- 
serves that the frequency becomes zero abruptly and the 
relaxation rate suddenly drops with a quadratic behavior 
when b c ~ 2kHz. The minimum of the parabola occurs 
at the magic angle, when 6 = 0. Then, all the polariza- 
tion reaching the 13 C at this orientation originates from 
protons outside the molecule. Then, the rate of 0.5kHz 
obtained at this minimum constitutes an experimental 
estimation of this mechanism. This has to be compared 
with the almost constant value of 1/(2tse) = 1/(2t0) ~ 
2.2kHz observed outside the magic angle neighborhood. 
This justifies neglecting the J-coupling and the direct re- 
laxation of the 13 C polarization through the dipolar in- 
teraction with protons outside the molecule. In the fol- 
lowing we describe our stroboscopic model that accounts 
for the "anomalous" experimental behavior. 



III. THEORETICAL DESCRIPTION 

A. The system 

Let us consider M coupled 1/2 spins with a 
Hamiltonian: 

n = n z + E + K + /2] , (3) 

i<j 

where Hz = Sj=i ^ ( w l + SiOi) If is the Zeeman energy, 
with a mean Larmor frequency cjl- The second term 
is the spin-spin interaction: 6y/o»j = is Ising, and 
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FIG. 1: (Color online) Spin swapping dynamics in 13 C- 1 H. (a): Experimental 13 C polarization in Fe(C5Hs)2 as a function 
of the contact time t and spin-spin coupling b(0). (b): Numerical simulations of the 13 C polarization obtained from Eq. J7J 
for different values of b, a dipolar system-environment interaction (\a/(3\ = 2) and a constant value for tse (tse = 0.275 ms) 
obtained by fitting the experimental data in the regime where the MKBE expression is valid. Projection plots in the b~t plane 
show a canyon where the oscillation period diverges indicating a Quantum Dynamical Phase Transition. 



ciij/bij = 0, 1, —2 gives an XY, an isotropic (Heisenberg) 
or the truncated dipolar (secular), respectively. This last 
case is typical in solid-state NMR experiments*^ where 
hSuJi, ctij,bij < hwh- 



In order to describe the experimental system let us 
take the first N = 2 spins, I\ (a 13 C) and / 2 (its di- 
rectly bonded 1 H), as the "system" where the swapping 
IJ.1T2) IT1J.2) occurs under the action of 612. The other 
M — N spins (all the other 1 H), with M — ► 00, are the 
spin-bath or "environment" . This limit enables the appli- 
cation of the Fermi Golden Rule or a more sophisticated 
procedure to obtain a mcanlife tse for the system levels. 
We will not need much detail for the parameters of the 
spin bath in Eq. J3J except for stating that it is charac- 
terized by an energy scale g?b which leads to a very short 
correlation time tb — h/ds- 



B. Spin^Fermion mapping 

The spin system can be mapped into a fermion parti- 
cle system using the Jordan- Wigner transformation) 2 ^^ 

if = cf exp |i7r^}=i c / c i}" Under the experimental 
conditions, Suji — 0, a 12 — and 6 i2 = b, the system 
Hamiltonian becomes: 

H s = fiwL {ctc x + c+c 2 - 1) + b (c+c 2 + c+ Cl ) /2. (4) 

The Jordan- Wigner transformation maps a linear many- 
body XY spin Hamiltonian into a system of non- 
interacting fermions. This leads us to solve a one-body 
problem, reducing the dimension of the Hilbert space 
from 2 N to N states that represent local excitationsi^Si 2 ^ 
To simplify the presentation, and without loss of gener- 
ality, we consider a single connection between the sys- 
tem and the spin bath, a 2 3 = ac? 2 3 and £> 2 3 = /W 2 3 with 
a 2j = b 2 j = 0, j = 4. . . 00 and oy = b±j = 0, j = 3 . . . 00. 



4 




blfi [KHz] 

FIG. 2: (Color online) Decoherence rate 1/(2t^) and fre- 
quency lo in the spin swapping of a 13 C- 1 H system. Data 
points are obtained by fitting cross polarization experiments 

MKBE 

to the expression P (t). The zero plateau in the fre- 
quency and the parabolic behavior of l/(2r^) in the region 
brsE <C h are indicative of an over-damped Zeno phase. Solid 
lines are the prediction of our model assuming a constant 
tse = 0.275 ms. 



Spins Ii with 3 < i < M are interacting among them. 
The SE interaction becomes: 

V = d 23 [a (c+c 2 - i) (c+c 3 -|)+/3 (c+c 3 + c+c 2 ) /2] . 

(5) 

In the first (Ising) term, the first factor "measures" if 
there is a particle at site 2 while the second "measures" at 
site 3. Hence, polarization at site 2 is "detected" by the 
environment. The hopping term swaps particles between 
bath and system. 

In the experimental initial condition, all spins are po- 
larized with the exception of IxM^ In the high tem- 
perature limit (hjJi,/kBT = s -C 1), the reduced den- 
sity operator is cr(0) = f G< (0) = (1 + s/f ) / Tr {1} 
which under the Jordan- Wigner transformation becomes 
Tr{i} ^ Tr{i} c ^ c 2- Since the first term does not con- 
tribute to the dynamics, we retain only the second term 
and normalize it to the occupation factor. This means 
that site 1 is empty while site 2 and sites at the particle 
reservoir are "full" . This describes the tiny excess above 
the mean occupation 1/2. To find the dynamics of the 
reduced density matrix of the "system" a (t) = ^G < (t) , 
we will take advantage of the particle representation and 
use an integral form^SiSi of the Keldysh formalism in- 
stead of the standard Liouville-von Neumann differential 
equation. There, any perturbation term is accounted to 
infinite order ensuring the proper density normalization. 
The interaction with the bath is local and, because of 
the fast dynamics in the bath, it can be taken as instan- 



taneous. Hence, the evolution is further simplified into 
an integral form of the Generalized Landauer Biittiker 
Equation (GLBE) 2 2*2I for the particle density. There, 
the environment plays the role of a local measurement 
apparatus. Firstly, we discuss the physics in this GLBE 
by developing a discrete time version where the measure- 
ments occur at a stroboscopic time T str .. 

C. Stroboscopic Decoherent Model 

We introduce our computational procedure opera- 
tionally for an Ising form ((3 /a — 0) of V. The initial 
state of the isolated two-spin system evolves with Hs- 
At T s tr. , the spin bath interacts instantaneously with the 
"system" interrupting it with a probability p. The ac- 
tual physical time for the SE interaction is then obtained 
as tse = T str . /p. Considering that the dynamical time 
scale of the bath (tb — h/d-o) is much faster than that 
of the system (fast fluctuation approximation), the dy- 
namics of site 3 produces an energy fluctuation on site 
2 that destroys the coherence of the two-spin "system" . 
This represents the "measurement" process that collapses 
the "system" state. At time T str . , the "system" evolution 
splits into three alternatives: with probability 1 — p the 
state survives the interruption and continues its undis- 
turbed evolution, while with probability p the system is 
effectively interrupted and its evolution starts again from 
each of the two eigenstates of cjc 2 . These three possible 
states at T str . evolve freely until the system is monitored 
again at time 2T str . and a new branching of alternatives 
is produced as represented in the scheme of Fig. 

A similar reasoning holds when (3 ^ 0. The sequence 
for isotropic interaction (a — (3) is shown in Fig. [3J}. The 
XY part of V can inject a particle. When an interrup- 
tion occurs, the bath "measures" at site 2 and, if found 
empty, it injects a particle. In the figure, this can be 
interpreted as a "pruning" of some incoherent branches 
increasing the global coherence. This explains why the 
rate of decoherence is greater when the Ising part of V 
dominates over the XY part. This occurs with the dipo- 
lar interaction and, in less degree, with the isotropic one. 
This contrasts with a pure XY interaction where the sur- 
vival of spin coherences is manifested by a "spin wave" 
behavior^ The empty site of the bath is refilled and 
de-correlates in a time much shorter than the swapping 
between sites 1 and 2 (fast fluctuation approximation). 
Consequently, the injection can only occur from the bath 
toward the system. 

D. Analytical solution 

The free evolution between interruptions is governed 
by the system's propagators Us (t) = exp[— iHst/h] = 
— yG 0R (i). The spin bath interacts with the system 
stroboscopically through an instantaneous interruption 
function (t) — ji- Is±£ -X < (t). The reduced density 
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FIG. 3: Quantum branching sequence for the swapping dy- 
namics. Panel (a) stands for an Ising system-environment 
interaction and (b) an isotropic one. Single states repre- 
sent states with interrupted evolution (incoherent) while pairs 
of states are coherent superpositions. Notice the self-similar 
structure. 



function is 



a(t)=V s (t)<r(0)U s (t)(l-p) n + 

n _ , 

E u s (t - t m ) s< (t m ) u f s (t - t m ) P (i - P y 



(6) 



where n = int (i/T s t r .) is the number of interruptions and 
t m = rar str .. The first term in the rhs is the coher- 
ent system evolution weighted by its survival probability 
(1 — p) n . This is the upper branch in Fig. [31 The second 
term is the incoherent evolution involving all the deco- 
herent branches. The m-th term in the sum represents 
the evolution that had its last interruption at mr s t r . and 
survived coherently from mr s t r . to nr str .. Each of these 
is composed by all the interrupted branches in Fig. |3| 
with a single state at the hierarchy level to with whom 
the paired state in the upper branch, generated at time 
(to + l)r s t r ., keeps coherence up to nr str .. 

In a continuous time process (rgs = T stI Jp with p — > 
and 7" s tr. — > 0), the above procedure gives a physical 
meaning for the Keldysh's self-energy 20,21 as an instan- 
taneous interruption function: S < (t) = Sjfj (t) + Sj^ (t) 
(see appendix). Dissipation processes are in Sj (t) while 



X< (i) involves only decoherence. The term (t) in- 
jects a particle, increasing the system density, provided 
that site 2 is empty. This occurs at an injection rate 
Pxy/tse = T XY /h with PXY = [P 2 / (a 2 + (3 2 )] , the XY 
interaction weight. The decoherent part, (t) , is a 
consequence of the "measurement" (or interruption) per- 
formed by the environment. It collapses the state van- 
ishing the non-diagonal terms (coherences) with a mea- 
surement rate Tzz/ft = (1 — Pxy) /tse- This also oc- 
curs with the XY interaction with rate Yxy/Ti regard- 
less of the fact that an effective injection takes place. 
Then, the inverse of the survival time or interruption 
rate is H/t se = Y SE = r zz + r XY . Unlike £< (t), 
this process conserves the total polarization. We denote 
by Y/h the rate associated to the isotropic interaction 
(a = [3 = 1). Then, the anisotropic rates are: Yzz = ct 2 Y 
and r XY = /3 2 r, then H/t se = Y SE = (a 2 + f3 2 ) Y. 
The rates can be calculated to infinite order in perturba- 
tion theory from the solution of the bath dynamics in a 
chain of spins with XY interactions^ which in the limit 
<W<fe -> gives r XY = 2tt (3 2 d 2 3 {nd B )- 1 = (3 2 Y which 
recovers a Fermi golden rule evaluation. 

The anisotropy ratio (a /ft) 2 accounts for the observed 
competition 3 ^ between the Ising and XY terms of V. 
The Ising interaction drives the "system" to the inter- 
nal quasi-equilibrium. In contrast, the XY term allows 
the thermal equilibration with the bath>2£ Solving Eq. 
© in the limit of continuos time processes we obtain 
our main analytical result for the swapping probability 
(experimental 13 C polarization): 



P(t) = l-a e 



-Rot 



- ai cos [(w + LRa) t + 4> ] 



-Rxt 



(7) 



which, in spite of appearance, has a single fitting param- 
eter. This is because the real functions u>, Ro, R\ and R2 
as well as ao, a\ and (j>o depend exclusively on b, tse and 
Pxy- Besides, b and pxy are determined from crystal- 
lography and the anisotropy of the magnetic interaction 
(pxy = 1/5 for dipolar) respectively. The phase transi- 
tion is ensured by the condition L0R2 = 0. The complete 
analytical expression is given in the appendix. 



E. Limiting cases 

Typical solutions of the quantum master equation 16 
for a spin swapping 30 follow that of MKBEi^. They 
considered an isotropic interaction with the spin envi- 
ronment, represented by a phenomenological relaxation 
rate 1/(2tse)- Within the fast fluctuation approxima- 
tion and neglecting non-secular terms, this leads to 

MKBE , . . - _ 

r (t) , used m most 01 the experimental fittings. 
Our Eq. (JJJ reproduces this result with 1/ (2r^) = 
R ~ 1/ (2tse) by considering an isotropic interac- 
tion Hamiltonian a = (3 = 1 under the condition 
1/ (2tse) <C b/h. However, at short times t <C tse, 
the MKBE swapping probability growths exponentially 
with a rate 1/ (2tse). In contrast, our solution manifests 
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the quadratic polarization growth in time, (^b/fi) t 2 . In 
Eq. this is made possible by the phase O in the 
cosine. In the opposite parametric region, 6tse -C h, 
our model enables the manifestation of the quantum 
Zeno effect ^S^S^ This means that the bath interrupts 
the system through measurements too frequently, freez- 
ing its evolution. At longer times, t 3> tse, one gets 
1 - P (t) ~ (l + | {b/hf r s 2 E ) exp [-\{b/KfT^t] , and 
the quantum Zeno effect is manifested in the reduction of 
the decay rate 1/ (2r<^) oc (b/fr) 2 tse as tse gets smaller 
than h/b. This surprising dependence deserves some in- 
terpretation. First, we notice that a strong interaction 
with the bath makes the 1 H spin to fluctuate, according 
to the Fermi golden rule, at a rate 1/t se . The effect on 
the 13 C is again estimated in a fast fluctuation approx- 
imation as 1/t0 cx (b/H) tse oc 

This "nesting" of two Fermi golden rule rates is formally 
obtained from a continued fractions evaluation of the 
self-energies 3 ^ 3 , involving an infinite order perturbation 
theory. Another relevant result is that the frequency de- 
pends not only on b, but also on tse- A remarkable differ- 
ence between the quantum master equation and our for- 
mulation concerns the final state. In the quantum master 
equation a (oo) must be hinted beforehand, while here 
it is reached dynamically from the interaction with the 
spin bath. Here, the reduced density, whose trace gives 
the system polarization, can fluctuate until it reaches its 
equilibrium value. 



F. Comparison with the experiments 



In order to see how well our model reproduces the ex- 
perimental behavior we plot the 13 C polarization with 
realistic parameters. Since the system is dominated by 
the dipolar SE interaction^ we take \a/(3\ — 2. We in- 
troduce b with its angular dependence according to Eq. 
(0 and we select a constant value for tse = 0.275 ms 
representative of the b 3> S/tse regime. Since in this 
work, we are only interested in the qualitative aspects of 
the critical behavior of the dynamics, there is no need to 
introduce tse(#) as a fitting parameter. These evolutions 
are normalized at the maximum contact time (3 ms) ex- 
perimentally acquired. They are shown in Fig. \T]p where 
the qualitative agreement with the experimental observa- 
tion of a canyon is evident. Notice that the experimen- 
tal canyon is less deep than the theoretical one. This is 
due to intermolecular 13 C- 1 !! couplings neglected in the 
model. We will show that the analytical expression of 
Eq. allows one to determine the edges of the canyon 
which are the critical points of what we will call a Quan- 
tum Dynamical Phase Transition. 



IV. QUANTUM DYNAMICAL PHASE 
TRANSITION 



Our quantum observable (the local spin polarization) is 
a binary random variable. The dynamics of its ensemble 
average (swapping probability), as described by Eq. J7J, 
depends parametrically on the "noisy" fluctuations of the 
environment through tse- Thus, following Horsthemke 
and Lefever^& one can identify the precise value for tse 
where a qualitative change in the functional form of this 
probability occurs as the critical point of a phase transi- 
tion. This is evidenced by the functional change (nonan- 
alyticity) of the dependence of the observables (e.g. the 
swapping frequency uj) on the control parameter br^/h. 
Since the control parameter switches among dynamical 
regimes we call this phenomenon a Quantum Dynamical 
Phase Transition. 

It should be remarked that the effect of other spins 
on the two spin system introduces non-commuting per- 
turbing operators (symmetry breaking perturbations) 
which produce non-linear dependences of the observ- 
ables. While this could account for the limiting dynami- 
cal regimes, it does not ensure a phase transition. A true 
phase transition needs a non-analyticity in these func- 
tions which is only enabled by taking the thermodynamic 
limit of an infinite number of spins^S In our formalism, 
this is incorporated through the imaginary part of the 
energy, S/tse, evaluated from the Fermi Golden Rulei 3 ^ 

When the SE interaction is anisotropic (a ^ (5), there 
is a functional dependence of u> on tse and b yielding a 
critical value for their product, 6tse/^ = fc PXY , where the 
dynamical regime switches. One identifies two paramet- 
ric regimes: I- The swapping phase, which is a form of 
sub-damped dynamics, when &tse/&p X y > ^ (^2 = in 
Eq. 0). 2- A Zeno phase, with an over-damped dynam- 
ics for &tse/&pxy < ^ as a consequence of the strong cou- 
pling with the environment (zero frequency, i.e. u) = 0, 
in Eq. 0). In the neighborhood of the critical point the 
swapping frequency takes the form (see appendix): 

bT S E/k pxY < h 

The parameters a PXY and k PXY only depend on pxy which 
is determined by the origin of the interaction Hamilto- 
nian. For typical interaction Hamiltonians the values of 
these parameters (pxy : k PXY , a PKY ) are: Ising (0, i,l), 

dipolar (|, 0.3564, 0.8755), isotropic (|,0,-^=) and XY 

(1,1,1). The swapping period is 




T, 



3/2 



(T 0c - T )- 1/2 . 



(9) 



where T = 2?2 is the isolated two spin period and its 
critical value Tq c = 2 , 7rrgE , determines the region where 

m 

the period T diverges as is observed in Fig. ^ The 
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estimated value of tse = 0.275 ms and dipolar SE inter- 
actions yield a critical value for the 13 C- 1 H coupling of 
b c /h = 2Tr/T 0c = 1.3 kHz. 

The complete phase diagram that accounts for the 
anisotropy of the SE interactions is shown in Fig. 0] 
There, the frequency dependence on pxv and brsv/h is 
displayed. At the critical line the frequency becomes zero 
setting the limits between both dynamical phases. 




FIG. 4: (Color online) Quantum dynamical phase diagram 
for the spin swapping operation. The figure shows the fre- 
quency dependence on system-environment (SE) interaction 
anisotropy pxy and the ratio among the internal and the SE 
interaction brsE/h. The projection over the forsB/ft — Pxy 
plane determines the phase diagram where the transition be- 
tween the swapping phase into the Zeno phase (uj = 0) is 
manifested. Values of pxy for typical SE interaction Hamil- 
tonian arc indicated in the contour plot. 



The two dynamical phases can now be identified in the 
NMR experiments which up to date defied explanation. 
The experimental setup provides a full scan of the param- 
eter 6tse through the phase transition that is manifested 
when the frequency goes suddenly to zero (Fig. |2t>) and 
the relaxation rate (Fig. changes its behavior decay- 
ing abruptly. The fact that 1/ (2t$) tends to zero when 
b -C S-/tse confirms the Zeno phase predicted by our 
model. In this regime, 1/ (2r^) is quadratic on b as pre- 
scribed. To make the comparison between the two panels 
of Fig. n quantitative, we fit the predicted dynamics of 
Fig. with P (f) . following the same procedure 
used to fit the experimental data. The solid line in Fig. 
|2]show the fitting parameters 1/ (2t^) and lu in excellent 



agreement with the experimental ones. 

We point out that Eq. @ is used to fit both the exper- 
iments and the theoretical prediction of Eq. 10 because 
it constitutes a simple, thought imperfect, way to extract 
the "outputs" (oscillation frequency u and a decoherence 
time r^). While the systematic errors shift the actual 
critical value of the control parameter, b/h, from 1.3 kHz 
to 2 kHz, Eq. @ yields a simplified way to "observe" 
the transition. 

V. CONCLUSIONS 

We found experimental evidence that environmen- 
tal interactions can drive a swapping gate through a 
Quantum Dynamical Phase Transition towards an over- 
damped or Zeno phase. The NMR experiments in spin 
swapping in a 13 C- 1 H system enable the identification 
and characterization of this phase transition as function 
of the ratio br^/h between the internal and SE interac- 
tion. We developed a model that describes both phases 
and the critical region with great detail, showing that 
it depends only on the nature of the interaction. In 
particular the phase transition does not occurs if the 
SE interaction is restricted to isotropic spin coupling. 
The phase transition is manifested not only in the ob- 
servable swapping frequency but also in the decoherence 
rate I/t^. While a standard Fermi golden rule pertur- 
bative estimation would tend to identify this rate with 
the SE interaction, i.e. = V t se — (<^23/^) 2 t b> 

as it occurs well inside the swapping phase, both rates 
differ substantially as the system enters the Zeno phase 
(fr^SE < k PXY h). Here the decoherence rate switches to 
the behavior cx (b/h) 2 t se . In the Zeno phase, the 

system's free evolution decays very fast with a rate t^ 1 . 
In spite of this, one can see that the initial state as a 
whole has a slow decay (its dynamics becomes almost 
frozen) because it is continuously fed by the environ- 
ment. Since the r SE has become the correlation time for 
the spin directly coupled to the environment, I/t^ pro- 
vided by our calculation can be interpreted as a "nested" 
Fermi golden rule rate emphasizing the non-perturbative 
nature of the result. Based on the wealth of this sim- 
ple swapping dynamics, we can foresee applications that 
range from tailoring the environments for a reduction of 
their decoherence on a given process to using the ob- 
served critical transition in frequency and decoherence 
rate as a tracer of the environment's nature. These ap- 
plications open new opportunities for both, the field of 
quantum information processing and the general physics 
and chemistry of open quantum systems. 3 ^ 
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APPENDIX A 

1. Discrete time process 

It is convenient to rewrite Eq. JfJJl within the Keldysh 
formalism where G < (t) = j-cr (t) and the evolution op- 
erator for the isolated system is G 0R (t) = — j-Us(i), 
where G 0A (t) — G 0R (t)^ , and the interruption function 
H < (t) ~ ■r^ I ^pX < (t) . We rearrange expression ©, for 

nT s tr. < t < (n + l)r s tr., in terms of nr str ., the last inter- 
ruption time: 

^G< (t) = G 0R (f - t n ) G< (t n ) G 0A (f - t n ) (1 - p) 



IG 0R (t - t n ) S< (t n ) G UA (t-t n )p 



OA , 



(Al) 



which takes advantage of the self similarity of the hierar- 
chy levels and is simple to iterate. It not only reproduces 
the results obtained from the quantum theory of con- 
tinuous irreversible processes22i2i discussed below, but it 
provides a very efficient algorithm which reduces memory 
storage and calculation time substantially. 



2. Continuous time process 

In the limit p — > and r str . — » with constant ratio, 
one gets a continuous process where t$e = T str ./p defines 
a survival time for the evolution of the isolated two-spin 
"system". Eq. (jHJl becomes 

G< (t) = h 2 G m (t) G< (0) G 0A (t) e~ t/TSE 

+ ( d^G 0R (t - t n ) £< (t n ) G 0A (t - t n ) e -(*-*")/^, 
Jo 

(A2) 

a new form of the GLBE22i2I that includes correlations 
(non-diagonal terms in G < and ^ < ). Moreover, keeping 
tse constant, under the conditions of finite T str . < fi/d.23 
and p < d,23/d,B (here oIb is a mean intra-bath interac- 
tion) there is no significant difference between the stro- 
boscopic approximation and the continuous time descrip- 
tion. Thus, the physics of a quantum theory of irre- 
versible processes^ is obtained by adapting a collapse 



theory^ into the probabilistic branching scheme repre- 
sented in Fig. The solutions of Eqs. 10 and (|A"2"|) 
are both computationally demanding since they involve a 
storage of all the previous time steps and reiterated sum- 
mations. However, Eq. I|AI|I provides a new computation 
procedure that only requires the information of a single 
previous step. Besides, it avoids resorting to the ran- 
dom averages required by models that include decoher- 
ence through stochastic or kicked-like perturbationsi^ZiiS 
Thus, it becomes a very practical tool to compute the 
dynamics in presence of decoherent and dissipation pro- 
cesses. 

The instantaneous interruption function: S < (t) = 
(0 + ^ M becomes, in matrix form, 



S< (t) 
S< (t) 



^G< (t) 





22 M 



■ >; /- ,"vy I 
Z1 TSE I [1 - f G< (t)] 



(A3) 
(A4) 



where the physical meaning of the Keldysh's self- 
energjiS2i2i becomes evident. 

Inserting Eq. (|A3|) and (IA4H in Eq. (|A2|I we obtain 
our main analytical result, Eq. Q, 



P(t) = iG< 1 (t) 

= 1 - a e~ Rot - a 1 cos [(<j + iR 2 ) t + <j) ] 



,-Rit 



(A5) 



here all parameters are real with 

= T~ r^IPxy,^) + 6— )b, (A6) 

2x \6 n(pxY,x)J 



where L0R2 = and are evaluated with x = oT^/h using 



4>{PXY- 



and 



V(PXY,x) = 

{4 (1 - PXY ) (to 2 - 2 (1 - pxy) 2 + 18p XY 
4-1 9 Is (at & - ((\ - pxY f + 12p XY \ ^ 



3 [4x 

4 Pxy ( 5 (1 - Pxy) 2 + 3p!y 



2 \ 2 

x 



- 4 Pxy ((! -Pxy) 2 -Pxy 



Also, 
Rq = 



6 



<ft(PXY, 
ViPXY, 



- Iv(pxy,x) + 
6 

I , 

Pxy + x (1 - pxy, , 

■"> / tse 



(A7) 



3 / 1 , A 1 
#1 = o Pxy + Pxy) 

I V / T SE 



i?0 

~2~ : 



(A8) 



and 



fl = 7T 



1 2 (w 2 - i? 2 ) + 2i? 2 - b 2 



2(a; 2 - J Rl) + (i? -ii i r' 
1 

02 = — ; :X 

2(uj + \R 2 ) 
(ggpgi - 6 2 ) (go - gi) + 2 (uj 2 - R 2 ) R 
(^-RD + iR.-R.f 

1 b 2 + 2R 2 - 4i? i? 1 

a-3 = r 



2(^_ jR 2 ) + (jRo _ jRi) 2' 
03' 



a 2 = a 2 , + a 2 , tan (0 O 



The parametric dependence of the swapping frequency 
uj, Eq. I|A6(1 . has a critical point at 2 = k PXY , where the 
transition occurs. In the neighborhood of this critical 
point uj and R2 takes the form: 



X Opxy) = -5291p XY - 1084p XY + 546p| Y - 4p XY + 1- 



24V3pxY V (28p XY - 2p XY + I) 3 , 



C (pxy) = -215p| Y - 2p XY + 1, 



and 



1 (A 2/3 + 36/ 2 ) (-Mi 2/3 + 36/ 2 / 3 + A/4 



A 5/3 A 



A = 36fc PXY - 8 + 24p XY + 48p| Y 



(All) 



- 36fc PxY pxY - 64p XY + 12/4, 
h = l/3fc PxY - l/3p XY " 1/9 (1 ~ Pxy) 2 , 



and 
where 



a 



, (b/h) - A: 2 XY /r g 2 E bT SE /k PKY > h 
bT S E/k PXY < h 



(A9) 



bT SE /k PXY > h 



ft,, 



pxyV fc lx Y / T SB ~ (V^) 2 br SE /k PXY < h ' 



fr 2 - 



1 

12 



(pxY - 1) x(pxy) ' +C(pxy) + 

- x 

(PXY - I) 3 C (PXY) 



19PXY 



[x(pxy)]' 



, (A10) 



A = - A + Pxy A - 6fc pxY - 16j4y + 20p XY 
— 10p XY + 13p XY fc PXY — 2k PKY pxY + k PXY , 



A = (l2fc PxY + 96p XY fc PxY - 120p XY fc PxY + 
^PxykpxY 39fcp XY p XY + 6fcp XY pxY 3fcp XY 

-48p XY + 48p XY -12p XY )^. 

In equation (|Allfl the functions Ai Aj A an( i A only 
depend of pxy- 
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